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Abstract 

We propose a new approach to the quantization of the damped harmonic oscillator in the 
framework of deformation quantization. The quantization is performed in the Schrodinger pic- 
ture by a star-product induced by a modified "Poisson bracket" . We determine the eigenstates 
in the damped regime and compute the transition probability between states of the undamped 
harmonic oscillator after the system was submitted to dissipation. 

1 Introduction 

Many attempts to study dissipative systems from the quantum point of view have been made for 
many years. Nevertheless, the problem in general is far from having a satisfactory solution. 

Quantum Mechanics mainly deals with microscopic systems and since dissipation is a macro- 
scopic concept some authors pointed out that the quantization problem is not even relevant because 
the nature of interactions or forces involved in such processes are not fundamental 1 . However 
there are problems where the dissipation at this level has an important role as in the quantum 
analysis of the radiation field in a microwave cavity, in quantum optics and also for deep inelastic 
processes in heavy ion scattering, to mention some of them ^ [2] . 

The damped harmonic oscillator (DHO) is one of the simplest cases of this type treated in the 
literature (for a review on the quantization problem see [3j and the references therein) . Most of the 
attempts have been using the canonical quantization formalism, dealing either with Hamiltonians 
dependent or independent of time, or by duplicating the number of degrees of freedom to obtain a 
Hamiltonian system that can be quantized by standard methods. We believe that the quantization 
of the DHO from first principles should bring new insights and ideas for a better understanding of 
the quantum aspects of dissipative phenomena. 

In this work we propose a quantization of the DHO from the point of view of deformation 
quantization. Our approach consists to include all of the dissipative effects into the algebra of 
observables. This is achieved by defining a star-product that takes into account the dissipation 
character of the DHO. We follow an idea of Enz 0j to add a symmetric part to the Poisson bracket 
(mathematically, a Hochschild coboundary) and work with the star-product * 7 induced by this 
modified bracket. 
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2 E-mail addresses: fturrub@u-bourgogne.fr, fturrub@fis.cinvestav.mx 



Chruscihski [SJ |S] has also considered the quantization problem for dissipative systems within 
the deformation quantization framework (see also [7]). His treatment differs from ours in the sense 
that he deals with the Moyal product and we rather use a mathematically equivalent star-product. 
In the limit when the damping factor goes to 0, we recover the Moyal product. 

Our paper is organized as follows. In Section 2 we review fundamental notions of deforma- 
tion quantization and give some details on the quantization of the harmonic oscillator in both 
Schrodinger and Heisenberg pictures. We introduce in Section 3 a new bracket between functions 
by adding a "dissipative" symmetric part to the canonical Poisson bracket. The modified bracket 
induces a ^-product that provides a quantization of the damped harmonic oscillator and allows 
obtain the corresponding quantum states. In Section 4, we compute some relevant transition 
probablities and observe that the dissipation couples a given eigenstate to lower energy states. 

2 Deformation quantization and the harmonic oscillator 

In this section we briefly introduce the basic notions of the deformation quantization formalism. 
We review the quantization of the harmonic oscillator in this setting by putting the emphasis on 
the Moyal-Schrodinger picture that will be used later on for the damped case. 

2.1 Notions on deformation quantization 

Deformation quantization was introduced in [S| as an alternative approach to quantization. In this 
formalism, quantization is interpreted as a deformation of the algebra of smooth functions on the 
classical phase space. We refer the reader to the original papers jS] for details and the review [5] 
for more recent developments. 

The phase space of a classical system is given by a manifold M endowed with a Poisson bracket 
P on the algebra of complex- valued smooth functions C°°{M). We shall denote by C°°(M)[[?i]] 
the space of formal series in H with coefficients in C°°(M). An element / in C°°(M) [[7i]] can be 
written as: 

oo 

/ = !>*/* w 

k=0 

where the fk's are in C°°(M). Rather than the Planck constant H, we shall use for physical reasons 
ih/2 (i 2 = —I) as deformation parameter. 

Definition. A deformation quantization of a classical system defined by (M, P) is an associative 
algebra (C7°°(M)[[fi]],*) oser C[[^]], where the associative product *, the so-called star-product, 

oo 

f*g = Y J {iK/2) k C k {f,g), f,geC°°(M), (2) 

fc=0 

satisfies: 

a) Co(f,g) = fg, and the Ck 's, k > 1, are bidifferential operators on C°°(M), 

b) the constant function I is the unit for *, i.e., !*/ = /*! = /, for any smooth function f , 

c) C\{f,g) — C\{g, f) = 2P(f,g), so that the star-bracket [/,<?]* = (f*g — g*f)/iti is a deformation 
of the classical Poisson bracket P. 

It is a famous result of Kontsevich ^U] that any Poisson manifold (M, P) admits a deformation 
quantization. 

The algebra of classical observables O is the set of real-valued functions on the phase space. 
Given a classical Hamiltonian H S O, the quantum evolution of an observable / is still governed 
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by the classical Hamiltonian, but one replaces the Poisson bracket by the star-bracket in order to 
get the Heisenberg equation of motion: 

|- = [/, HI = f * H 7 h H * f = P(f, H) + corrections in h. (3) 

As h — > 0, one recovers in some sense the classical equation of motion. The heuristic solution to 
Eq. © is given by f t = Exp^ (=$f) * f * Ex Pi (!§) , where 

v in ' ' n\ in' 

n=0 

with H* n = H • • * H (n factors). The star-exponential J3J is in general a quite singular object, 
nevertheless, for many physically relevant examples, the series converges as a distribution (see 
[S]). An important feature of the star-exponential is that its Fourier transform (when it exists in 
a distribution sense): J e~ ltE l h dv(E) allows to identify the spectrum of H as the support of the 
measure v. More precisely, suppose for simplicity that the Hamiltonian is such that the series in 
Q defines a periodic distribution in t, then the measure v is atomic and the Fourier transform 
reduces to a Fourier series (i.e. discrete spectrum): 

ExP ^ll^ =f>^„. (5) 

n— 

The -En's are the eigenvalues of H and the functions 7r„ represent the orthonormal eigenstates in 
the sense that: 

H * 7r„ = E n n n , n m -kTr n — 5 mn n n . (6) 

Strictly speaking this makes sense if one is using a convergent form of the product *. One may 
think of 7r„ as the function sent by the quantization rule to the projection operator |n)(n| in 
the Heisenberg picture of Quantum Mechanics, {\n)} being the normalized eigenstates of the 
Hamiltonian operator H corresponding to H. 



2.2 The harmonic oscillator in the Moyal-Heisenberg picture 

Let us consider the Hamiltonian H(q,p) = -^p 2 + 1J1 ^-q 2 of the harmonic oscillator on the phase 
space R 2 , with coordinates (q,p), endowed with the canonical Poisson bracket 

ptfg\^^L^l-^L^l (7) 
dq dp dp dq 

The Moyal product between two functions f,g £ C°°(R 2 ) is given by the formal series: 

f*g = exp (f P)(/, 9 ) ee fg + £ I(|) fc p fc (/, 5 ), (8) 

k>l 

where Pk is the & th -power of P, defined by 

P k (f,9Mp) ^ ^^)* /(ft,fll)fl( *' Pa) |n^- (9) 

The star-exponential for the harmonic oscillator has been computed in |S] and the series @ 
converges as a distribution to 

Exp*(-) = (cos (T )) ex P (— tan( y) ). (10) 
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The Fourier expansion of the right-hand side of i|l(J|) is given by [5j : 



tH 

Exp^(— ) =^exp(-i(n + l/2)u;t)7r n , (11) 

71=0 

where 

■K n (q,p) = 2exp ( — )(~l) n L n { — — — ), (12) 

and where L n is the Laguerre polynomial of degree n. In particular, the Gaussian function 

M q ,p)=2e X p C 2H n ^ p) ) (13) 

describes the vacuum state and, as expected, the energy levels of H are E n = huj(n+ 1/2). All this 
is well known and gives the essence of the deformation quantization of the harmonic oscillator in 
the Heisenberg picture. We shall below consider the less-known Schrodinger picture in deformation 
quantization. 

2.3 The harmonic oscillator in the Moyal-Schrodinger picture 

In order to write the Schrodinger equation: 

in dm =H + m (14) 

it is implicitly assumed that we have picked up some representation of the algebra of observables 
(0,*), so that we can 'multiply' observables with states. We refer the reader to the excellent 
review by Waldmann and references therein for a comprehensive exposition on states and 
representations in deformation quantization. 

We shall consider quantum states as a class of functions on the phase space. Strictly speaking 
this make sense only if we have a convergent (integral) form of the Moyal product * between a 
certain class of functions and h should no longer be a formal parameter. Indeed, solutions of 
the eigenvalue equation H * ip = Etp are not formal series in h, and it does not make sense to 
multiply such functions with the product ★ as formal series. Though we do not want to present 
a mathematical analysis of the convergence of the product ★ here, let us only mention that a 
convergent form does exist and gives a meaning to the product in such cases. 

For convenience, let us introduce the (dimensionless) annihilation and creation functions: 



, s P . muj p . muj , . 

aM = V^~l^ a{q ' p)= V^ + l^n q - (15) 

It is easily checked that we have: 

H = fkuaa = hw{a-k a + 1/2), [a, H]+ = — iua, [a, H]+ — iuia, [a, a]* = — . (16) 

The vacuum p in the Moyal-Schrodinger picture is the same as the one provided by the star- 
exponential i-e., po — ttoj an d in view of the relations Ijltjfl . it is clear that j5 n = a*- ■ - *a* p 
(there are n a's) solve the eigenvalue problem: H ★ p n — E n p n with E n = hw(n + 1/2). A simple 
computation using the Moyal product gives p n = 2 n a n po. 

In order to get a Hilbert space of states, we shall normalize the p n 's with respect to the scalar 
product: 

{f\9)tf = / f(<l,P)g(<liP) d K<liP)i (17) 
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where dp(q,p) = is the Liouville measure. A computation using Gaussian integrals shows 
that: 

Pn = -= a n Po , n > 0, (18) 
Vn\ 

is an orthonormal set for the scalar product (|17fl . 

The vector space over C generated by {p n }n>o is Sm = {h(a)TTo\h polynomial} and by taking 
the closure of Sm with respect to the scalar product (fTTjl . one gets the Hilbert space Hm = 
{f(a)TTo\f entire function such that J \f (a(q, p))\ 2 e~s^~ dp,(q, p) < oo}. Since H = huiaa, it is clear 
that TLm is isometric to the Bargmann-Fock space TLbf of antiholomorphic functions / on C such 
that: 

</I/>bf = - I \m\ 2 e-^m < oo. (19) 
T Jc 

Notice that Sm can be seen as a dense invariant domain for H in the sense that H ★ Sm C Sm- 
Any state ip can be expanded on the basis l|18|) : ip(q,p) = 53™>o a nPn{<liP), where a„ € C and 
Sn>o l a n| 2 = 1 m ^ uu analogy with standard Quantum Mechanics. 

3 The damped harmonic oscillator 

We shall introduce a new star-product * 7 that will provide a quantization in the Schrodinger 
picture of the damped harmonic oscillator whose dynamics is governed by: 

q + 2-/q + LU 2 q = 0, (20) 

where 7 > denotes the damping factor and uj is the frequency. 

Our philosophy is to incorporate all of the dissipative effects in the algebra of observables C 7 by 
an appropriate deformation of the Moyal product, with 7 as a second deformation parameter. We 
would like not to deform the measurable quantities and put everything related with dissipation into 
the star-product * 7 . For example, the function H = ^p 2 + 7J1 ^-q 2 still represents the mechanical 
energy of the particle submitted to dissipation. At the quantum level, we would like to keep H 
as the generator of some 7-deformed dynamics so that l|2U|l will still be deduced from the usual 
undamped Hamiltonian in a sense that we will now make precise. 

We mention that a singular Hamiltonian for the underdamped case has been derived in |12j . 
However the Hamiltonian considered there is only locally defined on the phase space. Consequently, 
the solutions of the induced system of equations are defined for small times only and cannot be 
extended globally in time. 

3.1 The modified bracket 

Consider M 2 as the phase space with coordinates (q,p). Eq l|2(J|l is equivalent to the system of 
equations: 

q = — , p — —muj 2 q — 2jp. (21) 
m 

We shall rewrite these equations in a form similar to the Hamilton equations of motion by using 
a modified bracket M obtained by adding a symmetric part to the canonical Poisson bracket P. 
This idea to consider a modified bracket for the quantization of dissipative systems within the 
Heisenberg picture is due to Enz 0|. 
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It is easy to check that the only bracket M that is a derivation in each of its argument and such 
that Eq. i|21|) can be cast into a form mimicking the Hamilton equations of motion: 

q = M{q,H) = £, 
m 

p = M (p, H) = -muj 2 q - 27p, (22) 

is given by 

M f,g Qq dp dp dq ^ m dp dp' 

Here is another motivation for this modified bracket and justification for the form of Eq. I|22l) . 
Let us consider a solution t i— > (x(t),y(t)) of Eq. (|21|l satisfying the initial conditions x(Q) = q 
and y(0) = p. The evolution of any smooth function / 6 C°°(R 2 ) is given by the path t i— > /t in 
C°°(R 2 ) defined by ft(q,p) = f{x(t),y(t)). The evolution equation of / t is deduced by taking its 
time derivative and using Eq. (|21H : 

lf(x(t),y(t)) = ^(x(t),y(t))i(t) + ^-(x(t),y(t))y(t) 

= ^ x (t),y(t))—-^-(x(t),y(t)) (m^x(t) + 2 ly (t)) 7 
oq m op 

= P(f,H)(x(t),y(t))-^rn^(x(t),y(t)) ^L(x(t),y(t)) 

= M(f,H)(x(t),y(t)). 

Thus the introduction of the modified bracket M appears to be quite natural. In particular, the 
time evolution of the mechanical energy H is 

H = M(H, H) = — —p(t) 2 < 0, (24) 
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and exhibits energy loss due to dissipation. 



3.2 The star-product * 7 

From the mathematical and deformation theory points of view, the bracket M has the following 
important property: it is a constant coefficient bidifferential operator of order 1 in each of its 
argument. This implies that M is a constant coefficient Hochschild 2-cocycle (see e.g. [5] for 
details on the Hochschild cohomology). As an immediate consequence of the previous fact, the 
following formula: 

f* 7 g = eMjM)(f,g) (25) 

defines a formal associative product 8 , i.e., a star-product. The formula (|25l) should be interpreted 

in the usual way: / * 7 g = fg + J2k>i h(t) Mk{f,g), where M k is the fc th -power of M defined 
in a similar fashion as in 

The star-product appears as a 7-deformation of the Moyal product *m m the sense that in the 
limit 7 — > + we indeed recover the Moyal product. Actually the relation between the product * 7 
and the Moyal product can be made explicit. First notice that the bracket M can also be written 
as: 

M(f, g) = P(f, g) + 2m^59(f, g), (26) 

where 8(f) = and 69 is the Hochschild differential of 8 defined by 69(f,g) = f9(g) — 

9(fg) + 9(f)g. It implies that the product * 7 is mathematically equivalent [5] to the Moyal 
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product, i.e., there exists a formal series of differential operators T = Id + X^iXt)^ 1, so ^ na * 
T(f *m 5) = T(/) * 7 T(g). Here the form J2g) of M gives us that: 

T = expHSm 7 e) = exp(^^^). (27) 

We shall use the product * 7 to quantize the damped harmonic oscillator. Notice that * 7 is not 
a Hermitian product f *-y 9 7^ 5*7 /j where the bar denotes complex conjugation. In fact 

the operator T destroys the Hermitian nature of the Moyal product as T(f) 7^ T(f) and this is 
precisely what is needed for dissipative phenomena. 



3.3 The quantum states 

The Schrodinger equation in our context reads: 



ih^-=H* 7 ^(t). (28) 

The remark after l)14|l also applies to the product * 7 . 

Hcuristically, the solution of the Schrodinger equation is given by: 

iP(t) = U(t) * 7 V(0), (29) 

where U(t) — Exp^(^-) is the star-exponential of H for the product * 7 . The function U(t) can 
be in principle determined by solving the differential equation: 

ih^- = H* 7 U(t), U(0) = 1, 

but it is easier to deduced it from the star-exponential i|10|) of H for the Moyal product and the 
use of the equivalence operator (|27[1 . A long but straightforward computation gives: 

AH. exp(7i/2) . , 

Exp, ( — ) = *y t x 30 

y *-> y ih> cosM/2)(l + ^tanM/2)) 1 /2 v > 

x exp ( — - tan(wt/2)(mo) 2 (7 2 H ^ )\ 

\huj K 1 JK q m(l + ^tan(wi/ 2 )) / 

According to the general principle of deformation quantization [Sj, the spectrum and the states for 
the Heisenbcrg picture are obtained by performing a Fourier decomposition of the star-exponential. 
Hence the periodic term in t in Eq. l|3Ufl . can be expanded as a Fourier series: 

Ex P ^(^)( g) p)=e Xp ( 7 t/2)^e XP (^)<( g ,p). (31) 

Again using the equivalence operator J27J, one finds that A„ = hw(n + 1/2) and thus the energy 
levels are given by = hu}(n + 1/2 + ij/2cu). Notice that the spectrum is now complex as 
a consequence of the fact that we are dealing with a nonhcrmitian product or, equivalently, an 
eventual operator representation of the algebra of observables for the product * 7 would represent 
H as a non-selfadjoint operator. 

The eigenstates TT^(q,p) in the Heisenberg picture can also be computed with the help of J23), 
but we just need to do it for the vacuum 7rJ, and it is given by: 

2 -1 p 2 
Ko{&P) = Ti o / u/2 e x P(-r-(Wy + —r, . , J )- (32) 
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(All over the paper, the branch cut of the square root is along the negative real axis.) As 7 — > + , 
we recover the vacuum 11130 for the Moyal case. It turns out that 7Tq is normalized for the I/ 2 -scalar 
product i|17|) : J |7Tq | 2 cZ/j,(g, j>) = 1 and it is annihilated by a(q,p), i.e., a * 7 7Tq = 0. Moreover we 
have: 

if = ftw(a* 1 a + 1/2 + i-y/2w), [a, H]* y = —ioja, [a, H]* y = icoa, [a, = (33) 
As in the Moyal case, it follows from the preceding relations that the fonctions: 

fL = a * 7 • • • * 7 a * 7 pi, {pl = ttJ) 

solve the ^-eigenvalue problem for H, i.e., H * 7 p% = fku(n + 1/2 + i-y/2u))p%. 

The /5^'s can be explicitly computed by taking the n th derivative in the variable s of the gener- 
ating function G(s,q,p) = J2n>o fn"Pra(s>P)- G(s,q,p) satisfies the differential equation: 



dG 
~d7 



i h dG 1 1 — - — / i2 7 ,5G 

aG —\ —Vmhiull H )— — , 

2V2 V mu dq 2V2 u) ' op 



with the initial condition G(0, q,p) = pQ(q,p), the solution of which is found by elementary means 
to be: 

G(s, q, P ) = exp (2sa(q, 3-^-) - ^^-j-)^). (34) 

The exponential factor in 1)34(1 . after a suitable scaling of the variable s, is the generating function 
of the Hermite polynomials H n : exp(2sz — s 2 ) = ^2n>o TI\Hn{z), and one easily find that: 

l — l2"f/u) l"f/LU l — l2j/U! 

The term of degree n in H n (z) is 2 n z n , and in the limit 7 — > + , we indeed recover p n = 2™a"poi 
the unnormalized states obtained before for the Moyal product. The normalized eigenstates for 
the scalar product lfT7|l are given by: 

pn/2 1 p 

where T = '7/% < and thus we have: 

* 7 ^ = fa;(n + 1/2 + * 7 /2^K, = S jf , (37) 

The Hilbert space 7i 7 for the damped harmonic oscillator is obtained in a natural way by taking 
the closure with respect to the L 2 -scalar product of the vector space generated by {p2,}- Let 
us introduce the complex- valued function b(q,p) whose complex conjugate is defined by b(q,p) = 
a(q, 1 _ i 2 1 /u )■ Hence 7i 7 consists of functions /(&(<?, J>))Po (9, p), where / is a holomorphic function 
in the variable b so that (/(&)Pq |/(^)Po)l 2 < 00 ■ Notice that H has a complex spectrum so it does 
not lead a self-adjoint operator on 7i 7 , one can check that (H -k (/>\ip) £2 7^ ((f)\H * ip) t 2 in general. 
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4 Measurements 



In the previous section, we have obtained an algebraic description of the dissipative effects of the 
damped harmonic oscillator by means of the product *~. In order to get a better understanding 
of the * 7 quantization during dissipation, we shall consider the following simple toy-model. 

Consider the usual one-dimensional harmonic oscillator. At time t — 0, one switches on the 
dissipation (7 > 0) and let the system evolves until time t = t > 0, at which time the dissipation 
is switched off so that we recover again the undamped case (7 = 0). 

Before t = 0, we have the usual harmonic oscillator described in terms of the Moyal star product 
*Af an d the eigenstates p n (q,p) given in ljT%|) . For < t < r, the quantum dynamics is governed 
by the product * 7 and the eigenstates p^(q,p) appearing in l|3tj[l . and for t > r the dynamics is 
given again by the Moyal product. 

Suppose that the system is known to be in the eigenstate p n (q,p) of the harmonic oscillator for 
t < (the phase is chosen so that it equals 1 for t = 0). This state belongs to the Hilbert space 
H.M defined at the end of Sect. 2.3. As the damping factor 7 is turned on at t = 0, the state of 
the system would then be in the Hilbert space TL 1 of the DHO. At t = 0, we use the following 
prescription 

Pn (q,p) = -$=(a(q,p)) n p (q,p) - -±=(b(q,p)) n p^p) (38) 



to express states in Hm as states in 7i 7 when the dissipation appears. Notice that the prescrip- 
tion H38J) sends the vacuum in Hm to the vacuum in 7i 7 , and boils down to the identity map when 
7 = 0. The vacua define the measure in the L 2 -scalar product, and l|38|) amounts essentially to 
make the substitution p — > p 7 = 1 _ i 2 J / LU f° r the momentum. 

Therefore, according to our previous prescription, the state in TC 7 that would describe the system 
at t = is: 

W{q,P,0) = -^=(b(q,p)) n p2(q,p). (39) 



It is not an eigenstate of H for the product * 7 , but it is an element of 7i 7 and therefore can 
be expanded on the basis p^ given by (|36|l in 7i 7 . Using the following expression involving the 
Hermite polynomials: 

I [n/2] , - 

2 n 2s k \(n-2k)V ' 

where [n/2] is the integer part of n/2, it is readily seen that the expansion of ip 7 (q,p 7 0) has the 
form: 

[n/2] 

r(q,P,0)=2Z a{ k n) Pl-2 k (l,P), (40) 

fc=0 



We obtain after a simple computation 



(n) 



2 n k\ V (n- 2k)V 



fc = 0,l,...,[n/2]. (41) 



Between t = and t = t, the state ip 7 evolves according to the dynamics defined by the product 



* 7 and we have: 



1/2] 



r(q,P,t) = £ 4 n) ex P (^-f^)p^ 2fc (g,p), < t < r, (42) 



ih 

k=0 
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where = ftuj(m + 1/2 + ij/2uj). 

When t — t the dissipation is turned off and the system will evolve as a harmonic oscillator in 
the Hilbert space %. We apply the inverse of the prescription |(3HJ|, i.e., 

-±=(b(q,p)) n pUq,P) - -±=(a(q,p)rp (q,p) (43) 



so that the state ip 1 at t = r becomes an element ip of the Hilbert space TCm . After applying the 
prescription (|43|l and using the eigenstates JSBJl, we find 

ln/2] M El -2k T Y n l 2 - k 1 
ip{q,P,T) = V Qfcexp( " 2fc )-—=== H n _ 2k (-— a{q,p))p {q lP ). (44) 

fe= o ih V("- 2fc ) ! 1 

Notice again that ip(q,p,T) is not an eigenstate of the Hamiltonian H for the Moyal product. It 
can be expressed in terms of the basis <|18|) of Hm with the help of the explicit expression of the 
Hermite polynomials 

Hn(x) = ^(-1) % •,,, (230"-^. 

A straightforward computation gives 

[n/2] 

V(g,P,r) = J2 Pt ] ^)pn-2k{q,P), (45) 

where 



P { k\ T ) = j^iy („ "L)! ex P(~^(" + l/2 + »7/2^)r)(exp(W) - l)\ (46) 

for fc = 0, 1, . . . , [n/2]. Notice that due to the nonunitary evolution during the damping regime, 
the state ip(q,p,r) is no longer normalized. In order to maintain the probabilistic interpretation, 
i/j should be normalized with respect to l|17|) . We denote by VP the normalized state: 

[n/2] [n/2] 

*(?,P,r) = — - J] ^ n) (r)Pn-2 & (g,p) s J] # n) (r)/>»-2*(«,p), (47) 

71 ' k=0 k=0 

where AT„(r) = (ELI 7 !? 1 ( r )| 2 ) 1/2 - For * > r ' the evolution of the state # is given as usual by: 

[n/2] 

*(?,p,t) = J2 P { k n) (r)eM-Mn-2k + l/2)(t-T))p n ^ k {q, P ), t > r. (48) 

In view of (|47|l . we notice that for generic values of r, the state of energy Hui(n + 1/2) we 
started with can decay to states of lower energy level n — 2k of same parity as n with transition 
probability |/3[™^(t)| 2 . Hence the state of energy E n spreads over lower energy states and the 
expectation value for the mechanical energy H would be less than E„ after the system has expe- 
rienced dissipation. However, for specific values of r, i.e. when lot = ir + 2nnr, m G N, we have 
/3^(t) = <5ofc exp(— iuj{n + l/2)r) and the system is still in the state p n and its behavior seems 
not be affected by the dissipation. 
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5 Conluding remarks 



Taking as a model the quantum description of the harmonic oscillator in deformation quantization, 
we have obtained a quantization for the damped harmonic oscillator in terms of a new star product 
•kn/. This product can be viewed as a deformation in two parameters, H and 7, that allows to recover 
at any point the classical case when h — > and the usual quantized harmonic oscillator when 7^0. 

The simple toy-model considered in Sect. 4 can be replaced by a more physical one. Instead of 
switching on/off abruptly the dissipation, a more physical situation would correspond to switching 
on/off the dissipation by letting 7 depend on time. 7(f) would then be a smooth positive function 
of t with support in some time interval [to,ti]. The star-product * 7 (t) would also depend on time 
and is defined exactly as in as the time-dependence of 7 does not alter the associativity. For 
a time s not in the support of j(t), * 7 ( s ) is simply the Moyal product and the whole system can 
be quantized with a unique time-dependent star-product. 

It is expected that the results obtained here will give more physical applications of the defor- 
mation quantization formalism and motivate the study of more complex systems. 
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